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Abstract 

We consider the Euclidean path integral approach to higher-derivative theo- 
ries proposed by Hawking and Hertog ( Phys. Rev. D65 (2002) 103515). The 
Pais-Uhlenbeck oscillator is studied in some detail. The operator algebra is re- 
constructed and the structure of the space of states revealed. It is shown that the 
quantum theory results from quantizing the classical complex dynamics in which 
the original dynamics is consistently immersed. The field-theoretical counterpart 
of Pais-Uhlenbeck oscillator is also considered. 
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I Introduction 



Theories described by Lagrangians containing higher derivatives seem to be of some 
relevance to physics. An interesting example is provided by Einstein gravity supplied 
by the terms containing higher powers of curvature [1] . Other examples include string 
theory with extrinsic curvature |2],[3], anyons [3], field theory on noncommutative 
space-time [4] or theory of radiation reaction [5] . 

The question whether a viable quantum counterpart of higher-derivative theory 
exists provides a real challenge. The straightforward canonical quantization based on 
Ostrogradski Hamiltonian formalism [6] suffers from serious drawback: the energy is 
unbounded from below. If one tries to manage this problem the ghost states seem 
to enter inevitably the theory; this is thought to be a fatal flaw (cf. [7j) unless it can 
be cured due to the existence of specific symmetry restoring unitarity in physical sector. 

It is sometimes claimed that the energy does not need to be bounded from below. 
In fact one can construct models which are stable in spite of the fact that their energy 
is unbounded from below. However, one should keep in mind that our models provide 
usually some idealization which consists in neglecting various perturbative terms that 
could easily spoil stability. On the other hand, there exists at least one general and 
physically relevant example of class of theories which contain ghosts but the unitarity 
in physical subspace is preserved for any physically acceptable perturbation; these are 
gauge theories. 

In general, in order to deal properly with the quantization problem one has to have 
a clear understanding of what is meant by correct quantization of a given classical 
theory. There are two natural conditions to be fulfilled: 

(i) the resulting structure posses all relevant properties of quantum theory ( super- 
position principle, reality of admissible values of measurable quantities, probabilities 
ranging between and 1, etc.); 

(a) in the limit H the initial classical theory should be recovered; in particular, 
the expectation values of basic variables for properly selected states (coherent states) 
obey in this limit the classical dynamical equations. 

The latter condition implies that the quantum theory is not only defined by the 
equations of motion and commutation rules but also by the transformation properties 
of dynamical variables under hermitean ( with respect to the relevant scalar product ) 
conjugation. This fact should be properly recognized when considering the admissible 
quantization schemes. It is often stressed that the Hamiltonian is not only defined by 
the formal differential operator but also by the boundary conditions imposed; this is 
fairly obvious. However, it is equally obvious that the latter cannot be imposed at will. 
In fact, the choice of boundary conditions influences the spectral properties of other 
dynamical variables (in particular - coordinates) which should in turn be chosen in 
accordance with the properties of their classical counterparts. For example, imposing 
in the case of harmonic oscillator the integrability condition along imaginary axis one 
obtains purely discrete negative spectrum. However, this is not a proper way to quan- 
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tize the classical harmonic oscillator as it contradicts the reality of coordinate variable 
(which is inherent in definition of harmonic oscillator). 

The celebrated Pais-Uhlenbeck (PU) oscillator |8j provides an even more distinct 
example. In the oscillatory regime it can be written ( within Ostrogradski framework 
) as a difference of two harmonic oscillators. Imposing the purely imaginary boundary 
condition on the second oscillator one finds the energy of PU system positive. How- 
ever, the original coordinate is a linear combination of oscillators coordinates and hence 
takes in general arbitrary complex values. On the other hand, the original coordinate 
of classical PU oscillator is real. Therefore, the choice of boundary conditions for the 
Hamiltonian enforcing the expected properties of the energy spectrum does not refiect 
the properties of other dynamical variables. If this is the case we are quantizing the 
classical theory which differs from the one we have started with; changing the initial 
theory is not a proper remedy for the problems we encounter. 

One can pose the question to what extent is the quantization procedure flexible; 
for example, whether some kind of complexiflcation is admissible. For standard La- 
grangians the dynamical equations are of second order in time derivatives. Therefore, 
classically the state of the system is determined uniquely once we know qi{t) and qi(t) 
at one given moment t. Any quantity of physical relevance (usually related to some 
Noether charge) can be constructed out of qi{t) and gi(t); moreover, the reality of qi(t) 
and qi{t) for a given t implies the reality of qi{t) and all its derivatives for any t. Thus 
qi and qi should be viewed as basic dynamical variables which, at the quantum level, 
generate the irreducible set of observables. Now, the problem is that, on the one hand, 
qi{t) and qi{t) are related by time derivative operation while, on the other, they can 
attain arbitrary values at a flxed initial moment. The way out is well known: it is 
the Hamiltonian formalism where qi{t) is replaced by pi which is related to qi{t) only 
"on-shell". Generalizing, one can say that the flrst step toward the quantization is to 
replace q^ by another variable in such a way that: (a) it is related to only "on-shell"; 
(6) the real character of q^ and (ji is preserved even if the additional variable takes 
complex values. 

In the case of Lagrangians containing also second order time derivatives the basic 
variables determining uniquely the classical state are q, q, q and Q . Therefore, accord- 
ing to the same line of reasoning one should immerse the initial dynamical system into 
a "Hamiltonian" one in such a way that q, q, q, Q are replaced by new variables related 
to the former ones only "on-shell". The new dynamics may be complex provided q{t) 
remains real for all t . So there is some freedom in the choice of the extended dynamics; 
the consistency condition, apart from reahty of q{t) is that q(t) obeys the dynamical 
equation(s) following from the initial action principle. Such an extended dynamics 
provides the starting point for quantization. 

The Ostrogradski formalism is the standard way of constructing the extended dy- 
namics. All variables, including the auxiliary ones, are here real and the theory is a real 
Hamiltonian one. Straightforward quantization leads to the Hamiltonian unbounded 
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from below and the space of states which is a Hilbert space (no negative norm states). 
One can also consider complex extensions of original dynamics. Such a possibility has 
been suggested some time ago [9], ^0\. It was found [9j to lead to positive energy and 
indefinite metric in the space of states. 

On the other hand, Hawking and Hertog [11] proposed the quantization method for 
higher- derivative theories based on Euclidean path integral in Lagrangian formulation. 
Our aim here is to study the Hawking-Hertog method in some detail. We compute the 
exact propagator for Pais-Uhlenbeck model and reconstruct the underlying quantum- 
mechanical algebra. The structure of the space of states is revealed and the existence of 
ghosts is shown to follow from the structure of classical action. Hamiltonian appears 
to be nonhermitean with respect to the standard scalar product. In spite of that it 
posses, in the oscillatory regime, purely real positive point spectrum. Outside this 
regime the eigenvalues cease to be real (at the crossover point the Hamiltonian is even 
nondiagonalizable) but the path integral continues to represent the propagator; the very 
existence of propagator is related to the fact that the real parts of energy eigenvalues 
are positive. We consider also the field-theoretic counterpart of PU oscillator. 



II The Pais -Uhlenbeck model 

In the Hawking-Hertog approach [11] one starts with the Lagrangian 

Lit)^L{q{t),q{t),m) (1) 
and assumes that its Euclidean counterpart 

LE{t) = -L{-it) (2) 
gives rise to positive-definite action 

t 

Se = j drLEir) (3) 



According to the suggestion of Hawking and Hertog one expects the Euclidean path 
integral J Dqe~^^^'^^ to converge and to give rise to a well-defined Euclidean quantum 
theory. 

Our aim here is to analyse the real-time counterpart of the Hawking-Hertog construc- 
tion for the Pais-Uhlenbeck Lagrangian: 

L=-q- — - y , m > 0, a > (4) 

We consider mainly the oscillatory regime 2ma < 1. The Euclidean version of eq.Q 
reads 

1 Tfl^ C? 
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We are interested in the propagator: 



- fdrLEir) 

K{q,q.;qo,qo;t) = I Dqe » (6) 



where the integration goes over all paths obeying q{0) = q^, g(0) = %, q{t) = q, q(t) = 
q. K is computed by standard method: first we split g(r) into classical path q^ obeying 
Lagrange equation together with the above boundary conditions and the quantum 
contribution to be integrated over. The result reads 

K = e-^^[«-'] / Dqe-^^^^^ (7) 



where 



t 

^Elq] = j drqi^r) ^- 



1 cP m . , , 

2^ + ¥ + T^'^^") 



The integral on the RHS of eq. (JTj) can be expressed in terms of the determinant of the 
differential operator entering eq.®. The latter is computed either by time discretiza- 
tion or by the method described in [10], up to a normalization factor to be chosen to 
provide the appropriate Dirac delta function for K in the t — )■ limit. The final result 
reads 

K{q, q- qo, %; t) = —^^e-^-t^-l (9) 

^^^^ ^2af sh^^yi + 2mat) _ sh\^Vl - 2mat ) 
ml 1 + 2ma 1 — 2ma 

Our aim is now to construct the quantum theory (observables, the space of states, etc.) 
in such a way that K, as given by eq. (Q, is the propagator in the standard coordinate 
representation. 

Naively, we could expect the following formula to hold 

K{q, q- go, go; i) = ^n{q. g)*„(go, go)e-*^" (10) 

n 

where '^n {En) are the eigenf unctions (eigenvalues) of the physical Hamiltonian H. Eq. 
( ITOl) implies the following symmetry properties 



K{q, q; go, go; t) = K{qo, go; g, g; t) (11) 
However, K, as given by eq.(l9]), obeys 



K{q, q; go, go; t) = K{qo, -go; q, -q; t) (12) 
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(actually, everything is real here). To see this let us note that, together with any 
solution Qci, qdi.'T) = Qdit — r) is also a solution to classical equation of motion obeying 
the boundary conditions: qd{0) = q, qci{0) = -q, qdif) = qo, qdif) = -qo] moreover, 
Sslqd] = Ssiqd]- 

By comparing eqs.( !T0|) ^( [T2|) one can conclude that the "physical" scalar product 
differs from the standard one ( cf. Ref. [9] and see below). 

The propagator K given by eq. (Q, obeys some differential equation. To derive it 
we check first, using standard methods [13.J that Sslqd] obeys the Hamiltonian-Jacobi 
equation with the Ostrogradski Hamiltonian He derived from Le, 

He = Piqi + - - ^Qi (13) 

with qi = q, q2 = q- This implies, by virtue of the identity 

C 2 d'SElqd] 



C a2 Qq2 

(which can be verified from the explicit form of 5*^;) that K obeys 
dK f d 1 92 1 



(14) 



together with the boundary condition K{q, q; go, qo] 0) = S{q — qo)S{q — qo) 
In order to construct the quantum theory from our Euclidean path integral we assume, 
as it is the case for first-order theory, that K is the kernel of the evolution operator 
K = e~*^ in the standard coordinate representation. K obeys 

^ = -HK (16) 
dt ^ ' 

Eq.f llSp is then identified with the coordinate representation version of eq. ffT6|l . This 
allows us to write the Hamiltonian 

1 Tn?' 1 
H = iq2Pi + -^pI + -Y^l + (17) 

Unfortunately, our Hamiltonian is not hermitean, ^ H. In order to cure this we 
introduce the metric operator rj obeying 

mi^ = (-i)''^^9i, vPiV = (-i)*^Vi, V = v^, v'^ = 1; (is) 

rj is defined by eq. (fT8|) up to a sign which can be fixed by demanding r]\qi,q2) = \qi,q2)- 
The new ("physical") scalar product is now defined by 

$) = (^, ?]$) = (^|?]|$) (19) 

One demands the observables to be hermitean with respect to the new product, i.e. 

A = A* = r]A^r] (20) 
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In particular H* = H. One can show that in spite of being hermitean only in 
generalized sense, has real positive spectrum and complete set of eigenfunctions. This 
can be formally proven by showing that one can convert H with the help of similarity 
transformation, into the sum of two harmonic oscillators 

H-B + + + (21) 

with Ai,2 = 2^(\/l + 2ma; =F \/\ — 2ma). However B, are defined only formally 
(at least one of them is unbounded so the domains should be carefully specified and 
shown to contain the relevant eigenvectors of the harmonic oscillators). A more precise 
algebraic proof can be given which is omitted here. 
So we infer the existence of complete set of eigenvectors '^ni,n2 obeying 

H'^n,,n, = (^(ni + ^)Ai + (n2 + ^)A2^ n^n, (22) 

Due to hermicity condition H* = H, \l'ni,n2 are mutually orthogonal with respect to 
the new scalar product; moreover, the normalization can be chosen in such a way that 

Let us expand any vector $ in terms of basic vectors \E'ni,n2 

oo 

^ ani,n2^ni,n2 (24) 
ni,n2=0 

Then, by virtue of eq.(l23l). 

«ni,n2 = (-l)"^(^ni,n2,$) = ("1)"^ (^«i,n2,^$) (25) 

and 

oo 

$=5^ (-l)"M^ni,n2,r/<f)vi>„,,„, (26) 

"1,12=0 

One can easily compute the new scalar product in coordinate representation: 

($,^) = j dqidq2{^il\qi,q2){quq2\^) = J dqidq2^iqi, -q2)'^{qi, ^2); (27) 
Let us now define the kernel K{qi, q2] qoi, qo2', i) by the equation 

(<l>2,e~*-^$i) = j dqidq2dqoidqo2^2{qi,-q2)K{qi,q2;qoiqo2;t)^i{qoi,qo2) (28) 
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Then 



($2, e-*^<l>i) = ($2, 77e'*^<l>i) = {r]^2, e~*^$i) = (29) 

J2i-^r (^$2,^n„nje-*^"-"^ (^ni,n2,r7$l) = 



ni,n2 

or 



Kiqi,q2;qoi,qo2;t) = J] (-l)"2e *^"i'"2^ni,n2(gi, g2)^ni,n2(goi, -^02) (30) 



ni,n2 



Now, K given as above, obeys eq. (IT^ and the appropriate boundary condition 
K{q, q; go, 0) = S{q — qo)S{q — Therefore, we conclude that the Euchdean path 
integral represents the evolution operator for the theory defined by the Hamiltonian f|T71) 
and the scalar product f|T9l) . 

It is worth to notice that our conclusions concerning nonhermicity of the Hamiltonian 
and the stucture of scalar product agree with those contained in Ref. [H]. 

Let us discuss the above results from the point of view presented in the Introduction. 
To this end consider the classical dynamics generated by the Hamiltonian (fT7|) . It reads 

qi = iq2, q2 = \p2 (31) 
Pi = -m'^qi, p2 = -q2 - ipi 
and results in the following equation for the variable q = qi 

o?qi^^ + ^1 + m^qi = (32) 

which is the Lagrange equation we have started with. Moreover, by inspecting eq. (1311) 
we find that it is consistent to assume that qi and pi are real while q2 and p2 are purely 
imaginary. The variables q = qi, q = iq2, q = -^P2 and q = ^^2 + characterizing 
uniquely the classical state, are real. Thus eq. fl3T]) present the consistent immersion 
of the initial dynamics in the extended, complex one; only the auxiliary dynamical 
variables take complex values. 

At the quantum level the expectation values of qi and pi are real; for §2 (as well as ^2), 
on the contrary, one gets 



($,^2$) = (g2$,$) = {^Al^) = -($,^2$) (33) 

because ^2 = ~?2- 

Let us note that the propagator can be also formally obtained from the Hamiltonian 
form of path integral provided the Hamiltonian ( 1T7|) is used (cf. also Ref. [10] ): 

t 

e-*^ = / DqiDq2DpiDp2e^ (34) 
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Another consistency check is provided by taking the t — )■ oo hmit of the propagator 
([9]). By virtue of eq. fl30|) one finds 



K{qi,q2-.qi,q2) ~ e^^^^^+^^^^^l/o o(gi, g2)^o,o(gl, -g^) (35) 
which can be easily verified to hold true. 

Few remarks concerning the Heisenberg picture are in order. Heisenberg equa- 
tions take the same form as the classical Hamilton equation (IHT]) . They preserve the 
*-hermitean (or *-antihermitean) character of dynamical variables; on the other hand 
the +-hermicity of canonical variables is not preserved: in general ql{t) ^ qi{t),pl{t) ^ 
Piif) even if qf{Q) = qi{0),pl{0) = Pi{0). However, this poses no problem since the 
standard Hilbert space structure plays only auxiliary role in our considerations; we 
could have started from the indefinite metric space from the very beginning. 

Up to now we were considering the oscillatory regime 2ma < 1. The natural 
question arises whether one obtains a reasonable quantum theory for 2ma > 1. By a 
closer inspection of the way eq.dH]) has been derived we immediately conclude that it 
is valid also in the region 2ma > 1. However, the preexponential factor in eq.(|n]) gets 
modified: 

Ory / sh'^(Tr-\/T+2rnat) si-n? ( tt- \/2ma— It) \ r ^ -i 



f^ff) = y m \ l+2mo 2ma—l J ('36') 

^ ^ \ 2a^sh\^) -e for 2m = I 

Eq.(|9]) shows that for 2ma > 1 the exponential damping is now modulated by oscillat- 
ing factors. Comparing this with the (assumed) form of K eq. f l30p . one concludes that 
the energies cannot be purely real. This conclusion is confirmed by the direct analysis 
of spectral properties of the Hamiltonian f|T7|) . It can be shown that for 2ma > 1 it 
posses the purely point complex spectrum 

ni + 712 /— — . ni -712 Y , . 

^ni.na = — V I + 2ma + I — V2ma - 1; (37) 

Za Za 



note that -E'ni,n2 = -^ni,ri2 S'^d as a result of reality of eigenequation, '^ni,n2 = ^n2,ni; 
moreover, ^'ni,n2 form a (nonorthogonal) basis in the space of states; 
For 2m.a = 1 the exponential damping is modified by the powers of t. Again, this 
agrees with the spectral properties of H: for 2ma = 1 its spectrum consists of the 
eigenvalues En = and to any n there corresponds an {n + l)-dimensional invariant 
subspace where H takes the Jordan block form (cf. [H] for other examples of nondi- 
agonalizable Hamiltonians) . 

We conclude that the Euclidean path integral gives the proper expression for the 
kernel of e~*^ also for 2ma > 1. In particular, the reality of K follows directly from the 
properties of the eigenvalues and eigenf unctions under complex conjugation. However, 
the resulting theory posses complex energy spectrum. 



9 



Ill Field theory 

Let us now consider the field-theoretic analogy of PU oscillator. We start with the 
Lagrangian 

1 1 7TI ^ -I- 711 ^ 

C = ^<l>(n + ml)(D + m2)$ = ^(D*)' + ^ ^ ^ <|.D$ + (38) 

2 2^2 I I, J. 2 ^ 2 

-m^mg'I' + o.t., > 

In analogy with eqs. ( [3l]) we put 

$ = vl>^ $ = 2Vl>2 = vi>^ $ = ivi/2 (39) 

and define the extended Lagrangian 

£(^,^,^2,^2,A) = (40) 

1 • TTi^ -\- Tfl^ ' 1 

_(,^2 - A^i)2 + ' ^i(2^2 - A^i) + -m?m2$2 ^ ^^^^ _ ^^^^ 

The theory is singular due to the appearance of Lagrangian multipliers A. Let 



L = / d^xC (41) 



The canonical momenta read 

6L 

TT = ^ = 
SX 



P2 = = -^2 + i ( -A^i + ^{mj + ml)) 



5^2 \ 

so there are two sets of primary constraints 



71 ^ 

Pi - A ^ (43) 



while the Hamiltonian density reads 



1 i 1 

n = --pI - ipzA^i + -{ml + m2)p2^i + - -m\m\^\ 

+ ^(m2 + m^)^2^p7rA + (T(pi -A) (44) 
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yielding the canonical equations 



X = p 
n = a 

^1 = + (45) 



^2 = -P2 - ^A^i + ^{ml + ma)^! 

pi = iAp2 - -^("^1 + "1-2)^2 - ^("T-i + "T'2)^^i + mlml'i/i 
h = -ipi 

Eqs.^ and ^ yield 
a = 

p = iAp2 — 2^^'^ ~^ "^2)^2 — ^("^1 + rnlY'^i + miml'i/i (46) 

Eqs.f l43p are the only constraints; they are of second kind and allow to eliminate A and 
TT. The Dirac brackets for remaining variables take the standard form. 
It is now straightforward to construct the quantum counterpart of our theory. We skip 
standard details and write out the final result: 

^ r 

$i(a;) = J2 (fk {ai{k)e~'^^^ + a*(A;)e*") (47) 
i=i ^ 



where 



Redefining 



^^Ck) = y"'^ ,,,, (49) 

we arrive at the (almost) standard commutation rules 

Mk),bl{k')] = (-1)^+15,^,5(3) (fc - fc') (50) 

The space of states is spanned by the fourmomentum eigenvectors describing particles 
of two types created by the operators aj(fc). The subspace of the positive norm consists 
of states carrying an even number of particles of type "2". 
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Let us consider the energy-momentum tensor. It is not difficult to find its symmetriezed 
form 



^-g'^''mlml<l>' (51) 



One easily checks that, in spite of its form, T'^^ depends on time derivatives of $ up 
to third order only. In particular, expressing $ = \E'i, $, $ and $ in terms of canonical 
variables (with the help of eqs.f l45p ) we find that T^^ coincides, up to some spatial 
boundary terms, with the Hamiltonian density (jS]). In terms of $ and its derivatives 
it reads 

yoo ^ 1^2 _ ^-^ _ i(A$)2 + $ A $ + 

^(m? + m^)($ A $ - $2) - irn^m^^^ (52) 
Inserting here the general solution ( H7|) one gets 

J d^xT^^ = j d?k {klhl{k)bi(k) + klbl(k)h2{k)^ (53) 

The energy has positive eigenvalues and conserves the particle numbers. The latter 
holds true also for other generators of Poincare group. Therefore, the subspace of 
positive definite scalar product is invariant under the action of Poincare group. The 
theory is perfectly consistent. 



IV Conclusions 

Let us compare the Hawking-Hertog approach with the standard one based on Ostro- 
gradski formalism. In both cases the dynamics of the original variable q is immersed 
in a wider one. In the standard case all additional variables are kept real while within 
the approach inspired by the Euclidean path integral some of them attain complex 
values. This difference results in difference in the form of the Hamiltonian. The Os- 
trogradski Hamiltonian is hermitean but unbounded from below; on the other hand, 
the Hamiltonian following from Euclidean path integral is nonhermitean which implies 
the necessity of modifying the definition of scalar product. In the oscillatory regime 
2'ma < 1 the nonhermitean Hamiltonian has real spectrum which, in contrast to that 
of Ostrogradski one, is bounded from below. For 2ma > 1 the Ostrogradski Hamilto- 
nian continues to be hermitean; its spectrum, however, becomes continuous (and still 
unbounded from below). On the other hand the Hamiltonian fll7p has a purely point 
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spectrum even for 2ma > 1 (while it is not diagonalizable if 2ma = 1) but the eigen- 
values are no longer real (although the real parts are positive). Therefore, it cannot be 
considered as a viable quantum mechanical Hamiltonian; its defects cannot be cured 
by any change of metrics in the space of states. It must be, however, stressed that the 
Euclidean path integral still represents the operator e~*^. 

Having a detailed knowledge of the structure of the space of states one can specify 
the kind of perturbations which generate the unitary evolution in the subspacc corre- 
sponding to the eigenvalue 1 of metric operator rj. They should be represented by the 
operators hermitean both in the usual as well as generalized sense. Then the evolu- 
tion operator commutes with rj and leaves invariant the above subspace. We should, 
however, assume that the interaction Lagrangian does not contain second derivatives; 
if this is not the case the theory will be probably not well defined. 

Similar results hold true for the field-theoretical counterpart of PU oscillator. We 
conclude that Hawking- Hertog method provides, at least in some specific 
sonablc, consistent method of quantization for theories containing second derivatives 
quadratically. 

The necessary condition for the applicability of H-H method is the positive defi- 
niteness of Euclidean action. However, as we have seen, this by far not sufficient: the 
Euclidean action for PU oscillator continues to be positive definite for 2ma > 1 while 
the resulting Hamiltonian has complex eigenvalues or is not diagonalzable. 
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